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The nonlinear time evolution of the quantum fields is studied in the O(N) model for large 
N in a radiation dominated FRW universe, with a view towards the phenomenon of explosive 
particle production due to either spinodal instabilities or parametric amplification, i.e. preheating. 
Quantum backreaction effects due to the produced particles are included consistently within the 
, large N approximation. We find that preheating persists when the expansion is included, although 

the amount of particle production is reduced compared to the values found in Minkowski space. We 
also see that the behavior of the evolving zero mode is very different from that in Minkowski space, 
though the late time behavior in all cases is determined by a sum rule that implies the existence of 
O i 1 Goldstone bosons in the final state. 

m : 

1 I. INTRODUCTION 

m : 

A viable model of inflation must be able to reheat the Universe to high enough temperatures to regain the tested 
features of the standard Big-Bang model, such as nucleosynthesis. However, it has recently been noted jl]-|3| that the 
\ theory of reheating as originally conceived [Q] requires revision due to abundant particle production originating from 
C\| ■ the coupling of the zero mode of the inflaton held to other modes of this same field or modes of a different field (see 
i also refs. §-|). 

Preheating is an inherently non-equilibrium process and as such, any field theoretical treatment must take this into 
account. Preheating is also a non-perturbative process; if A is either the self-coupling of the given field or the coupling 
Q\ • of this field to another, then, in typical scenarios for which the initial energy in the zero mode is 0(1/A), the number 
of particles produced is typically 0(1/A). Thus any perturbative treatment will necessarily be incomplete. 

In the Minkowski space analysis of refs. ||[|, these constraints were dealt with by applying the so-called closed 
time path formalism (CTP) flTo| of quantum field theory, which describes the non-equilibrium time evolution of the 
density matrix describing the state of the field theory to both the Hartree approximation of a quartically self-coupled 
scalar theory as well as the large N approximation to the O(N) vector model pj]JT^ |. These analyses clearly displayed 
the dissipational effects of backreaction due to the produced particles as the zero mode evolved towards its final 
equilibrium state. They also showed that the zero mode can execute extremely interesting and unexpected behavior, 
depending upon the couplings as well as the initial conditions ||. 

In this work, we ask the question: does the inclusion of the Hubble expansion affect the preheating process? Some 
partial answers to this question were given in ref. |l3|]l4|| . 

We consider the O(N) model embedded in a fixed radiation dominated FRW background spacetime. Again, we 
include the effects of backreaction to leading order in 1/N and calculate the behavior of the zero mode, as well as 
the particle production per mode. Our results are that the expansion does not prevent an epoch of preheating from 
occurring, but will in general reduce the number of particles produced per physical volume compared to the flat 
spacetime results. Furthermore, the energy lost from the zero mode due to redshift will typically prevent some of the 
unexpected behavior found for the zero mode in Minkowski space |,|| from occurring in an expanding universe. 



II. THE FORMALISM AND THE MODEL 



The correct way to describe the evolution of quantum fields in non-equilibrium situations is the CTP formalism 
pof which describes the time evolution of the density matrix of field theory. Details of the applications to field theory 
are given in ref. ||||. 

A controllable yet non-perturbative approximation scheme is the large N approximation of the O(N) vector model, 
which has been generalized for non-equilibrium situations | pT| . The lowest order term in this approximation corre- 
sponds to mean field theory but, at least in principle, the next to leading terms that would include effects such as 
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scattering, could be computed. Such corrections will be important in any discussion of how the particles produced 
during preheating actually thermalize, and in particular for any estimates of the reheating temperature. 

The 0(N) model in an FRW universe can be formulated in the following manner |15|. The Lagrangian density is: 



C = a 3 (t) 



V(<7,tt) = -^m 2 (t 



l- 2 l (V$(£,f)) 2 a 

2^ (M) ~2 a{t) 2 ~ V ^- 

+ — r (6-6f 



(2.1) 



for fixed A in the large N limit and m 2 > (symmetry breaking). Here a(t) is the FRW scale factor, and 4> = (<r, 7?) 
is an O(N) vector where if represents the N — 1 "pions" . For simplicity, we will only consider the case of flat spatial 
sections, and we will set the renormalized coupling to the Ricci scalar R to zero. 
We decompose the a field as 



a{x,t) = a (t) + x(x,t). 



(2.2) 



The large N limit can be implemented consistently with the following Hartree-like factorization (for more details the 
reader is referred to ref. B): 



N-l 



n(x,t) =V(*,t) (1,1, •••,!) , a (t)=4>(t)VN, 



(2.3) 



with 



(V/)«0(1), <x 2 )«0(l), 0*0(1). 



(2.4) 



The leading contribution is obtained by neglecting the 0(1/ N) terms in the formal large N limit. 

If Vfe(i) are the mode functions (for comoving wavenumbers k) corresponding to i/)(x, t), then the evolution equations 
for <j>(t) and the Vk(i) are: 



dt 2 



cf>(t) + 3H<t>(t) + <p{t) 
d k 2 



3H 



dt a 2 (t) 



A 



A 



+ o^(*)+o(r(*)> 



= 



2 + ^ 2 (t) + ^(t)) 



V k {t)=0, 



(2.5) 
(2.6) 



W k = 



d 3 k \V k {tf 
(2tt) 3 2W k 



(2.7) 
(2.8) 



where H = a(t)/a(t) is the Hubble parameter. The initial conditions are Vfc(to) = 1 , Vk(to) 
After renormalization flla.01, we have: 



-iWk 



V"(r) + 3^- V '(t) + M 2 (r) v {r) = 0, 



a(r) 

d 2 | ,fl'(r) d | g 2 
dr 2 a(r) dr o 2 (t) 



+ M 2 (t) 



V q (r) = 0, 



(2.9) 
(2.10) 



M 2 (t) = -1 + V 2 (t)+ 9 E(t), 



(2.11) 



where r = mt, q = k/m, T) 2 (r) = X(f> 2 (t) / 2m 2 , and g = X/8tt 2 . Primes denote derivatives with respect to r. The 
fluctuations operator gfS(r), which accounts for the back reaction of the produced particles on the modes, is: 



5 s (t) = g dq q 



W a 



1 



0(q 



qa 2 (r) 



2q 3 



1 



M z {t) - -R(t) 



a' 2 (to) 
a 2 (r) 



(2.12) 



2 



where k is a renormalization scale, R(t) is the Ricci scalar in units of |m 2 | and V q (ro) = 1 , V^(t ) = —iW q . 

We consider: (i) slow roll initial conditions with 77(70) ~ 0. In this case some of the modes display spinodal 
instabilities which then drive the production of particles; (ii) chaotic initial conditions, T)(to) 3> 1. The particle 
production in this case is due to parametric amplification of the non-zero modes due to the evolution of the zero 
mode. In both these cases, we will take as our background a radiation dominated FRW universe, with scale factor 

1/2 

given by a(t) = {{t + t )/t ) v . This analysis can be repeated for a matter dominated universe with similar results, 
while the analysis of the de Sitter case is quite different |L6| . 

A. Slow-Roll Initial Conditions 

We take 77(10), to be very near the origin: 77(70) = 10~ 5 , rf{ro) = 0. The coupling is fixed at g = 1CP 12 as would 
befit an inflationary theory. Larger couplings will tend to speed up the dynamics, but will leave the qualitative aspects 
unchanged. 

Figure la Figure lb Figure lc 




FIG. 1. (a) The zero mode rj(r) vs. r, (b) the quantum fluctuation operator D(t) vs. t, and (c) the number of pions N q vs 
q at r = 30 (dashed line), r = 150 (solid line) for the parameter values rjo = 10~ 5 ,tj' = 0,(7 = 10 _12 ,/io = 0.1. 

Fig.(|l]a) shows the behavior of r)(r) for an initial (dimensionless) Hubble parameter ho = 0.1. We see that it makes 
one excursion out to near the classical turning point and then tries to return to its initial value. However, by this 
time, the expansion of the universe has redshifted enough energy from the zero mode that it cannot come back to the 
origin and begins to execute damped oscillations about the tree-level minimum at 77 = 1. 

Fig.(|l|b) shows the quantum fluctuations operator D(t) = ffS(r) during the evolution. We see that D(r) grows due 
to spinodal instabilities as 77(7") makes its first excursion; this is when particle production occurs. However, at late 
times D(t) tends towards zero. 

Fig.([j]c) shows the number of "pions" per unit physical volume, in units of 1/g, produced as a function of comoving 
dimensionless wavenumber q at various times. An important feature to notice here is that this distribution has support 
at q = Q. 

In Minkowski space where hg = 0, the zero mode finds a stable minimum near the origin, while the fluctuations 
grow to D(t — > 00) — > 1. Furthermore, the momentum distribution of the produced pions has essentially no support 
near q = in this situation. Thus the addition of expansion induces a qualitative change in the dynamics of the 
system and the power spectrum of fluctuations. 

Increasing ho yields essentially the same features, the only differences being in the time scale involved. Thus for 
larger ho, 77 begins to oscillate about the tree-level minima sooner, and the fluctuations begin their decrease towards 
zero sooner. 

B. Chaotic Initial Conditions 

We take 77(70) = 4 , 77' (to) = ,ho = 0.1 as our chaotic initial conditions, again with g = 10~ 12 . In this case, the 
zero mode has enough energy to be able to sample both minima for a period of time until enough energy is lost both 
due to the redshifting as well as to particle production to restrict it to small oscillations about one of the tree level 
minima, in this case the one at 77 = 1 (Fig.(||a)). This behaviour is different than the Minkowski case || where the 
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zero mode keeps oscillating over the two minima showing that the symmetry is unbroken if the initial energy is larger 
than V(r) = 0). 



Figure 2a 





Figure 2c 
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FIG. 2. (a) The zero mode 77(7") vs r, (b) the quantum fluctuation operator D(t) vs. t, and (c) the number of pions N q vs 
q at r = 40 (dashed line), r = 200 (solid line) for the parameter values rjo — 4, r)' — 0, g — 10~ 12 , ho = 0.1. 




As in the slow-roll case, the fluctuations are again driven to zero at late times as the zero mode reaches its minimum 
(Fig.(Sb)). However, the behavior of the momentum distribution of particles (Fig.(||c)) is quite different, looking much 
more like the Minkowski slow-roll space distribution than the slow-roll one above. 

Changing the initial value of 77(70) will change the number of times both ground states are sampled; for 77(70) > 4, 
we have more oscillations back and forth before enough energy is drained away to make the zero mode settle into one of 
the vacua. If the initial Hubble parameter is increased, for fixed 77(70), the energy of the zero mode is redshifted away 
faster, so that fewer large amplitude oscillations are executed and fewer particles per physical volume are produced. 



III. LATE TIME ANALYSIS 



An interesting by-product of our numerical analysis concerns the late time behavior of the field. The damping term 
in the zero mode equation suggests that there will be stationary solutions at late times. Such solutions must obey 
the large N Ward identity 



?7oo(-l + vlo +5^00) = 0, 



(3d) 



where the subscript 00 denotes the late time value of the relevant quantity. Our numerical results confirm that eq. ( pd| ) 
is satisfied to one part in 10 10 . The same sum rule was found to be satisfied in the Minkowski case as well. For ijoo 7^ 
the term inside the bracket must vanish, leading to a mode equation for massless Goldstone bosons, which can be 
solved exactly in a radiation dominated universe: 



V q (r) 



1 



a(r) 



c„ exp 



(2 i9 T 1/2 (7 



TO 



d/2 



1/2 

d q exp(-2igT ' (t + t 



,1/2 



(3.2) 



where c q and d q are determined by the initial conditions and the dynamical evolution. For q 7^ and late times, 
r ^> 1/((? 2 to), the phases average out so that |V g (t)| 2 ~ a(i")~ 2 and the fluctuation contribution is driven to 



C 



TO 



oc a(r) 2 , 



(3.3) 



where C is a constant which depends on the dynamics of the evolution and initial conditions. Using cq^Sd), the zero 
mode then behaves as: 



T -C 



TO 



(3.4) 



If 77(70) 0, we find that Vo(r) tends to a constant at late times. This explains the appearance of the peak in N q 
at q — 0. Since 77(70) = , 77' (tq) = is a fixed point of the zero mode equation, eq.( |3.lD is satisfied with gEoo = 1- In 
this case only the q — > modes contribute substantially to the fluctuations and the power spectrum of the fluctuations 
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becomes both peaked at zero momentum and narrower at longer times. This analysis reproduces our numerical results 
to high accuracy. 

There is a striking contrast between the radiation dominated and the de Sitter cases. In the latter, for t)(tq) =/= 0, 
there are consistent stationary solutions for which is not a minimum of the tree level potential and the fluctuation 
contribution gHoo is driven to a nonzero constant fll6fl. 



IV. CONCLUSIONS 



We have generalized our previous results on the quantum dynamics of the O(N) model in Minkowski space to FRW 
universes. In particular, we find that preheating can still occur in an expanding universe, though the actual amount 
of particle production depends sensitively on the size of the initial Hubble parameter relative to the mass of the field. 



We also showed that the late time dynamics is determined by the sum rule in eq.(3.1). This, together with the 
behavior of the modes, drives the system to satisfy the sum rule by damping the fluctuations to zero and allowing the 
zero mode to find one of the tree level minima. This indicates that the symmetry is broken at late times JB| . We see 
no evidence of symmetry restoration J^]. 

The next step in this program is clearly to attempt to allow the scale factor to evolve dynamically and understand 
how the back reaction of the produced particles influences the dynamics of a(t). We are currently setting up the 
formalism to tackle this problem. 
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